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Abstract 

We discuss the interplay between spectral shape and detector response beyond a simple E'"^ 
neutrino flux at neutrino telescopes, at the example of time-integrated point source searches 
using IceCube-40 data. We use a self-consistent model for the neutrino production, in which 
protons interact with synchrotron photons from co-accelerated electrons, and we fully take 
into account the relevant pion and kaon production modes, the flavor composition at the 
source, flavor mixing, and magnetic field effects on the secondaries (pions, muon, and kaons). 
Since some of the model parameters can be related to the Hillas parameters R (size of the 
acceleration region) and B (magnetic field), we relate the detector response to the Hillas 
plane. In order to compare the response to different spectral shapes, we use the energy 
flux density as a measure for the pion production efficiency times luminosity of the source. 
We demonstrate that IceCube has a very good reach in this quantity for AGN nuclei and 
jets for all source declinations, while the spectra of sources with strong magnetic fields are 
found outside the optimal reach. We also demonstrate where neutrinos from kaon decays 
and muon tracks from r decays can be relevant for the detector response. Finally, we point 
out the complementarity between IceCube and other experiments sensitive to high-energy 
neutrinos, at the example of 2004-2008 Earth-skimming neutrino data from Auger. We 
illustrate that Auger, in principle, is better sensitive to the parameter region in the Hillas 
plane from which the highest-energetic cosmic rays may be expected in this model. 
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1 Introduction 



Neutrino telescopes, such as IceCube [l] or ANTARES [2j , are designed to detect neutrinos 
from astrophysical sources. If protons are accelerated in these astrophysical objects, as we 
expect from the observation of the highest energetic cosmic rays, the collision with target 
photons or protons will lead to charged pion production, and therefore to an extraterrestrial 
neutrino flux. Therefore, neutrino telescopes are an indirect method to search for the origin 
of the cosmic rays. There are numerous source candidates, the most prominent extragalactic 
ones being Active Galactic Nuclei (AGNs) [3]-|6| and Gamma- Ray Bursts (GRBs) [7J, see 
Ref. [s] for a review and Ref. for the general theory of the astrophysical neutrino sources. 



In generic estimates or bounds, such as Ref. |10, 11 , the cosmic ray flux is related to 
the potentially expected neutrino flux. Very interestingly, IceCube-40, referring to the 40 
string configuration of IceCube, is currently starting to touch these generic estimates for 



particular source candidates, see Ref. 12 for GRBs. In addition, time-integrated 13] and 



time-dependent 14 point source searches have been performed, so far, without success. 



On the other hand, there may be sources for which the optical counterpart is absorbed. 



so-called "hidden sources", see, e.g., Ref. 15 18 . This immediately raises questions how 



to interpret the data (see also Ref. 19 for AGN models), in particular: what does it 
mean that IceCube has not seen anything? What parts of the parameter space is IceCube 
actually most sensitive to? In this study, we address these questions in terms of the interplay 
between spectral shape expected from the sources and the detector response. In order to 



quantify the detector response, we use the time-integrated point source analysis in Ref. 13 
for different source declinations. In addition, we compare the parameter space coverage 
to other experiments, at the example of 2004-2008 Earth-skimming neutrino data from 
Auger [20]. 



A convenient description of the parameter space of interest is the Hillas plot 23 . In order to 
confine a particle in a magnetic field at the source, the Larmor radius has to be smaller than 
the extension of the acceleration region R. This can be translated into the Hillas condition 
for the maximal particle energy 

£^max [GeV] ~ 0.03 ■7]-Z-R [km] ■ [G] . (1) 

Here Z is the charge (number of unit charges) of the accelerated particle, B is the magnetic 
field in Gauss, and rj can be interpreted as an efficiency factor or linked to the characteristic 
velocity of the scattering centers. Potential cosmic ray sources are then often shown in a 
plot as a function of R and B, as it is illustrated in Fig. [T]by the numbered disks (see end 
of figure caption for possible source correspondences). Assuming that a source produces the 
highest energetic cosmic rays with E ~ 10^° eV, one can interpret Eq. ([T]) as a necessary 
condition excluding the region below the dashed line in Fig. [l] (for protons with t] = 0.1). 
However, this method does not take into account energy loss mechanisms, which may lead 
to a qualitatively different picture - as we shall see later. In the following, we will study the 
complete parameter space covered by Fig. [l] without any prejudice. Since the location of 
the sources in Fig. [T] cannot be taken for granted, we will refer to the individual sources in 
Fig.[l|as "test points" (TP) in most cases, and leave the actual interpretation to the reader. 

Our main focus is the particle physics perspective, i.e., we start off with the minimal ingre- 
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Figure 1 : Maximal proton energy as a function of R and B in our model (contours) for an acceleration 



efficiency rj = 0.1 on the Hillas plot; data taken from Fig. 2 in Ref. 21 . The dashed line indicates the Hillas 
condition in Eq. ([l]) for 10^° eV protons. The dotted line separates the regions where synchrotron (above 
line) and adiabatic (below line) energy losses dominate the maximal proton energy in the model, i.e., the 
maximal proton energy follows the Hillas condition below the dotted line. The region "UHECR" indicates 
where 10^° eV cosmic ray protons are expected to be produced in the model. The "no acceleration" region 
refers to either inefficient acceleration, or inefficient pion production. The different "test points" (numbered 
disks) correspond to: (1) neutron stars, (2) white dwarfs, (3) active galaxies: nuclei, (4) active galaxies: jets, 
(5) active galaxies: hot-spots, (6) active galaxies: lobes, (7) coUiding galaxies, (8) clusters, (9) galactic disk, 
(10) galactic halo, (11) supernova remnants, (12) additional test point. Test points taken from Ref. [22] . 

dients to neutrino production and the main impact factors to the spectral shape. First of 
all, note that most neutrino telescope analyses (and even some models) use an E'"^ neutrino 
spectrum as an initial assumption, see, e.g., Ref. |l3], which is often believed to be consistent 
with a proton injection spectrum (oc E'"^) coming from Fermi shock acceleration. However, 
in p7 interactions, the neutrino spectrum follows the pion spectrum, which depends on pro- 
ton and photon spectral shape. Therefore, the assumption for the neutrinos only holds 
for an E~^ target photon density (for the E'"^ proton injection spectrum), as it is often 
assumed for the prompt emission from GRBs up to the photon break. In self-consistent 
models, where the target photons originate, for example, from synchrotron emission of co- 
accelerated electrons, such a hard spectrum is very difficult to obtain, which results in a 
different neutrino spectrum; see, e.g., Ref. for an example. In this work, we study the 
detector response to different shapes of the neutrino spectra. For that purpose, we require a 
toy model which can predict the neutrino spectral shape for wide regions of the parameter 
space without any bias from a multi-messenger observation, since, after all, neutrino sources 
may not be necessarily seen as photon sources. The model needs to take into account the 
dominant particle physics processes especially affecting the spectral shape, which are multi- 



pion processes in the pion production (see, e.g., Refs. 24-26]), neutron and kaon production 
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see, e.g., Refs. 21,27-29]), magnetic field effects on the secondary muons, pions, and kaons 



see, e.g., Refs. 21,25,29-32]), and the appropriate maximal energies. Note that for the 



neutrino spectral shape, especially the magnetic field is an important control parameter, 
see, e.g., Refs. 21 33 for quantitative discussions. First of all, the maximal proton energy, 
which is recovered in the neutrino spectrum, can be limited by proton synchrotron emission. 
Second, for high enough magnetic fields, both the maximal energy of the neutrino spectrum 
and the spectral shape will be determined by the synchrotron cooling of the secondaries. 
Therefore, in order to predict the neutrino spectral shape, it is, especially for large B, more 
important to accurately model the particle physics of the secondaries instead of the cooling 
and escape processes of the primaries (protons, electrons, photons). Especially if the target 
photons come from synchrotron emission, there is little sensitivity to the spectral shape of 
the parents (only the square root of the parent spectral index enters the photon spectrum). 



We use the model in Ref. 21 for the prediction of the spectral shapes, where one of the 
starting points was the minimal set of assumptions required to describe the neutrino flavor 
ratios and spectral shapes as accurate as possible, given the above boundary conditions. 
In this model, charged pions are produced from photohadronic (^7) interactions between 
protons and the synchrotron photons from co-accelerated electrons (positrons). The photo- 
hadronic interactions are computed using the method described in Ref. |26j, based on the 
physics of SOPHIA [34], including higher resonances, t-channel charged pion production, 
and multi-pion production. The helicity-dependent muon decays are taken into account 
as described in Ref. |25]. The toy model relies on relatively few astrophysical parameters, 
the most important ones being the size of the acceleration region {R), the magnetic field 
strength at the source (B), and the injection index (a) which is assumed to be universal for 
protons and electrons/positrons. Naturally, the parameters R and B can be directly related 
to the Hillas plot. The leading kaon production mode and the energy losses of all secon- 
daries (muons, pions, kaons) are taken into account. Our key argument goes as follows: for 
each set of R, B, and a, we can predict the spectral flux shape using this model. Using the 
exposures from IceCube or Auger, we can then compute the sensitivity limit (normalization 
of the flux) for exactly this parameter set. In order to compare sensitivities to different 
spectra, the (integrated) neutrino energy flux density at the detector [ergcm~^s~^] is used. 
We illustrate at the end of Sec. [3| that this quantity is typically a direct measure for the 
luminosity x pion production efficiency of the source. Thus, "good sensitivity" in this work 
means that the source will be found for relatively small values of this product. Note that 
we do not interpret the data in terms of the deeper astrophysics involved, and we do not 
predict the absolute levels of the neutrino fluxes. In addition, note that our conclusions 
will naturally be somewhat model-dependent. However, as we will illustrate, some of the 
qualitative results should be rather robust. 

This study is organized as follows: In Sec. |2j we summarize the key ingredients of the 
model, where a critical discussion of the limitations of the model can be found at the end 
of this section. Then in Sec. |3j we describe our method to limit individual fluxes using the 
neutrino effective areas from IceCube-40. Furthermore, in Sec. |4]we show the result for the 
whole Hillas plane and different source declinations, we illustrate the impact of the injection 
index a, neutrinos from kaon decays, and z/T--induced muon tracks. In Sec. [5| we emphasize 
the complementarity of different experiments and data sets, using Auger as an additional 
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example. Finally, we summarize in Sec. [6j 



2 Review of the source model 



The model in Ref. |2l] describes neutrino production via photohadronic (p7) processes for 
transparent sources (optically thin to neutrons) and includes magnetic field effects on the 
secondary particles (pions, muons, kaons). It can be used to generate neutrino fluxes as a 
function of few astrophysical parameters. Below we outline the key ingredients of the model 



relevant for this study; for details see Refs. 21 , 26 , or Ref. 35 for a shorter summary. All 



of the following quantities refer to the frame where the target photon fleld is isotropic, such 
as the shock rest frame (SRF). 

The protons and electrons/positrons are injected with spectra cx E~"', where a is one of the 
main model parameters. The maximal energies of these spectra are determined by balancing 
the energy loss and acceleration timescale given by 

t-^\ = V— , (2) 

with 7] an acceleration efficiency depending on the acceleration mechanism, where we typi- 
cally choose ?7 = 0.1. If synchrotron losses dominate, the maximal energy is therefore given 
by 

^max - Y ■ 

It scales oc m^, which means that the protons are accelerated to much higher energies, and 
oc l/\fB, which means that strong magnetic fields limit the maximal energies. If adiabatic 
energy losses dominate, t~l ~ c/i?, the maximal energy is (for protons) given by Eq. ([l|, 
i.e., the Hillas condition. Since the neutrino energy follows the proton energy, the maximal 
energy of the protons determines (for not too strong magnetic field effects) the peak of the 
neutrino spectrum in E'^dN/dE. We show the maximal proton energy in Fig. [l]as a function 
of R and B for this model, where the dotted line separates the regions where synchrotron 
(above line) and adiabatic (below line) energy losses dominate the maximal proton energy 
(here rj = 0.1)]^ One can easily see that in the upper part (above the dotted line) the Hillas 
condition, Eq. ([T]), does not apply, since the synchrotron losses dominate, whereas below 
the dotted line, the maximal proton energy follows Eq. g. One can also get an educated 
guess for the best IceCube spectral shape sensitivity already, since its differential limits are 
minimal at around 10^ to 10^ GeV (see Sec. [s]) - which should be a very robust prediction. 
In the absence of magnetic field effects, one can estimate that the optimal detector response 
is roughly obtained when these differential limit minima coincide with a certain fraction of 



^Note that these assumptions are consistent with Ref. 36 (efheient acceleration case) and Ref. 37 , who 



have emphasized these effects earher. EspeciaUy, the dotted hne approximately corresponds to the middle 



solid curve in Fig. 6 of Ref. 37 , which limits the region where the Hillas condition for the maximal proton 
energy applies. That, however, does not mean that the region on the r.h.s. of the dotted line is excluded, 
it only means that the maximal energy is limited otherwise (by synchrotron losses, and, in the very lower 
right corner of our plot, also by interactions with the CMB). Compared to Ref. 37 , we show in Fig. [l] the 
maximal energy directly. 
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the maximal proton energy in Fig. [T] (the neutrinos take about 0.05 to 0.1 of the proton 
energy in photohadronic interactions). Note that in this model, the potential sources of the 
highest-energetic cosmic ray protons are to be found in the lower right corner in the region 
labeled "UHECR" within the 10^^ GeV contour (unless there is a strong Lorentz boost of 
the source). In this case, the neutrino energies extend up to about 10^° GeV. However, 



recent results on the cosmic ray composition, such as from Auger 38 , indicate that, at 
least above about 3 10^ GeV, heavier nuclei may dominate the cosmic ray composition. 
Since these heavier nuclei may reach the same energies for lower magnetic fields according 
to Eq. ([T]), they would occupy a slightly different region in the Hillas plane. Therefore, 
one should keep in mind that the mapping into the Hillas plane only exactly applies to the 
proton contribution at ultra-high energies. 

For each particle species, the injection and energy losses/escape are balanced by the steady 
state equation 

Q{E) = ^{h{E)N{E)) + ^, (4) 

with tesc(-E) the characteristic escape time, h{E) = -Et^-^l^ with t'l^iE) = -1/EdE/dt the 
rate characterizing energy losses, Q{E) the particle injection rate [(GeVscm^) ^] and N(E) 
the steady particle spectrum [(GeVcm'^) ]. For all charged particles, synchrotron energy 
losses and adiabatic cooling are taken into account. In addition, unstable secondaries, i.e., 
pions, muons, and kaons, may escape via decay. As a consequence, for pions, muons, and 
kaons, neglecting the adiabatic cooling, the (steady state) spectrum is loss-steepened above 
the energy 

where synchrotron cooling and decay rates are equal. One can read off this formula that the 
different secondaries, which have different masses m and rest frame lifetimes tq, will exhibit 
different break energies Ec oc \/m^/To which solely depend on particle physics properties 
and the value of B. These different break energies will lead to a spectral split of the neutrino 
spectra, which is an imprint of the magnetic field. 

While being accelerated, the electrons loose energy into synchrotron photons, which serve 
as the target photon field. Charged meson production then occurs via 

P + 7 TT + p' , (6) 

p + 7 -> ir+ + A/E , (7) 

with these synchrotron photons, where the leading kaon production mode is included and p' 
is a proton or neutron. In addition, two- and multi-pion production processes are included 
(not listed here), see Ref. I26j for details. The injection of the charged mesons is computed 



from the steady state proton Np{E.p) and photon N^{b) spectra with Ref. 26 

oo oo 

QbiEt,) = j ^Np{Ep) j deN,{e)Rh{x,y), 
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with X = Eh/Ep the fraction of energy going into the secondary, y = {Epe)/mp (directly 
related to the center of mass energy) and a "response function" Rb{x, y) (see Ref. ^^). The 
weak decays of the secondary mesons, such as 

7r+ ^ l^^ + iy,,, (9) 
/i+ ^ e+ + z/e + I^;, , (10) 

are described in Ref. [25] , including the helicity dependence of the muon decays. These will 
finally lead to neutrino fluxes from pion, muon, kaon, and neutron decays. 

In order to compute the z/^ (/3 = e, /i, r) neutrino flux at the detector, i.e., including flavor 
mixing, we sum over all these initial neutrino fluxes of flavor weighted by the usual flavor 
mixing 

3 

-Pa/3 = l^m.n^ffil^ , (11) 

1=1 

where Uai are the entries of the Pontecorvo-Maki-Nakagawa-Sakata mixing matrix. We use 



sin^ 9i2 = 0.318, sin^ 623 = 0.5, and sin^ 613 = for the sake of simphcity (see, e.g., Ref. 39 ), 
leading to flavor equipartition between z/^ and at the detector. We sum over neutrinos 
and antineutrinos, i.e., if we refer to "z/^", we mean the sum of the z/^ and z>^ fluxes. 

Concerning the limitations of the model, it certainly does not apply exactly to all types of 
sources. For example, in supernova remnants, pp (proton-proton) or pA (proton-nucleus) 
interactions may dominate the neutrino production, which would require additional param- 
eters to describe the target protons or nucleons. In addition, at ultra-high energies, heavier 
nuclei may be accelerated, see discussion above. The spirit of this model is different: It 
is developed as the simplest (minimal) possibility including non-trivial magnetic field and 
flavor effects. The relevant point for this study is that the instrument response depends 
on the neutrino energies and spectral features. The spectral features, as a peculiarity of 
neutrinos, are in cases of strong magnetic fields dominated by the cooling and decay of the 
secondaries (pions, muons, kaons). The energy of the secondaries follows the proton energy, 
to a first approximation, Ei, ~ x Ep with x = (9(0.1) to 0{1). These two observations 
apply to any of these mentioned interactions, no matter if pj, pp, or pA, if the energy of 
the accelerated particle is much larger than the one of the interaction partner (Feynman 
scaling). Therefore, we expect qualitatively similar results for these cases. 

Another variable is the target photon density, which is assumed to come from synchrotron 
emission of co-accelerated electrons here. In more realistic models, typically a combination 
of different radiation processes is at work. For example, there may be thermal photons 
radiated off an accretion disk, proton synchrotron photons, inverse Compton up-scattered 
photons, photons from pair annihilation, and photons from 7r° decays cascaded down to 
lower energies. It is therefore a frequently used approach for a particular type of source to 
compute the neutrino spectrum from a target photon spectrum which corresponds to the 
observation without describing the origin of these photons, see, e.g., Ref. [12j for gamma- 
ray bursts, where a band function parameterization of the observed gamma-ray fluxes is 
used. Qualitatively, the target photons control the interaction rate of the primaries in pj 
interactions [see, e.g., |26j], which enters the normalization and is not used here. In addition. 



6 



the spectral shape depends on the target photon spectrum in the energy range relevant for 
the interactions. In many examples with strong magnetic fields a spectral break in the 
photon spectrum is less important than the cooling and decay of the secondaries, which 
depend on particle physics only. For example, in our model, different hypotheses for the 
cooling of the electrons, leading to the target photons, have been tested, and in most cases 
the relevant part of the spectral shape hardly changes. Thus, while it is unlikely that the 
model applies exactly to a particular source, it may be used as a good starting hypothesis. In 
addition, note that from Eq. ([S]), only the product of the proton and photon (and therefore 
electron) density normalizations enters the final result. The final normalization of the 
neutrino spectrum will depend on the source luminosity, the interaction volume, a possible 
Lorentz boost of the acceleration region F, and the redshift z of the source. Since we only 
need the spectral shape, we do not compute the normalization explicitely. Note, however, 
that also the final neutrino energies depend on F and z. For the sake of simplicity, we neglect 
these effects, but one should keep in mind that the actually observed energy spectrum could 
be significantly shifted in energy as oc F/(l + z). 

The main parameters of the model are i?, affecting the shape and maximal energy (via 
^ad ) of primaries, 5, affecting the maximal energy (via t^nchr) primaries and the 

break (via of the secondaries, and a, affecting the spectral slope of the primaries. The 
parameters R and B can be directly related to the Hillas parameters, see Fig. [TJ 



3 Method and impact of the spectral shape 



Here we describe our method, i.e., how we constrain individual fluxes from IceCube data 
and quantify the response of the instrument. The simplest possible approach is total event 
rate-based, which means that no information on the reconstructed neutrino or muon energy 
is used explicitely. It can be described with the exposure Exp(i?, 5) = A^^{E^ 5) tcxp, where 
A^^ is the neutrino effective area and texp is the observation time. Here A^^[E^6^ is a 
function of the flavor or interaction type (which we do not show explicitely), the incident 
neutrino energy and the declination of the source 6. The neutrino effective area already 
includes Earth attenuation effects and event selection cuts to reduce the backgrounds, which 
depend on the type of source considered, the declination, and the assumptions for the input 
neutrino flux, such as the spectral shape. Normally, the cuts are optimized for an E~'^ 
flux, which means that for speciflc fluxes with different shapes, the following analysis may 
slightly improve by an optimization of the detector response. On the other hand, one has 
to understand that the experiments cannot optimize their event selection for any possible 
input spectrum. The total event rate of a neutrino telescope can be obtained by folding the 
input neutrino flux with the exposure as 

/V = / ,lE^ME,i) = / iEAf[E,i)uJ-^ . (12) 

Here dN{E)/dE is, for point sources, given in units of GeV~^ cm~^ s~^. If backgrounds are 
negligible, the 90% (Feldman-Cousins) sensitivity limit K^q for an arbitrarily normalized 



input flux used in Eq. (12) can be estimated as K^^q ~ 2AA/N [40]. This imples that 



a predicted flux at the level of the sensitivity limit, irrespective of the spectral shape. 
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would lead to the same number (2.44) of events. The 90% confidence level differential limit 
E'^dN/dE can be defined as 2.3 E/Exp{E, 6), see, e.g., Ref. We will comment on the 

impact of the backgrounds below. 

In the following, we use the neutrino effective areas for time-integrated point source searches 
in IceCube-40 [13], based on 2008-2009 data with texp = 375.5 days. We show these neutrino 
effective areas for muon tracks and for six declination bands, as in the original reference, in 
Fig. [2] For the sake of readability, we use "downgoing", "upgoing", or "(quasi-)horizontal" 
as additional description of the declination bands. Note that for IceCube located at the 
South Pole, a particular source will always appear under a specific declination, which means 
that the position in the sky will have some impact on the sensitivity, and the declination 
bands correspond to different source classes or data sets. Since upgoing muon neutrinos 
become absorbed in Earth matter above about 30 PeV and strong cuts have to be applied 
for the downgoing events to reduce the cosmic muon background, quasi-horizontal events 
imply the optimal performance. 

The effective areas shown in Fig. [2] are for muon tracks only. However, note that there are 
two contributions: not only produce muon tracks, but also Ur, since the tau lepton decays 
into a muon with a branching ratio of about 17%. We show the individual contributions 
in Fig. [2| where thin solid curves represent the z/^ effective areas, the dashed curves the z/^ 
effective areas, and the thick solid curves the total effective areas for either i/^ or u^- if flavor 
equipartition v^:Vt-=1:1 is assumed. Note that the total effective areas for the individual 
or fluxes are given by the maximum of the individual effective areas, not the sum. 
Assuming flavor equipartition, one can read off from Fig. |2] that quasi-horizontal muon 
tracks are dominated by z/^ events. The downgoing muon tracks also come mostly from z/^, 
but there is a small region at high energies where Vj. dominate. In this case, the muon tracks 
are recovered at lower energies. The upgoing muon tracks above 30 PeV, where the Earth 
becomes opaque to z/^, are dominated by z/^ events. Therefore, it is a priori not clear if and 
when the muon tracks from Vr can be neglected. In general, our analysis will be based on 



z/^ interactions, but we will come back to this point in Sec. 4.4 Note that beyond the limits 
shown in Fig. |2| we use linear extrapolation, which turns out to be a good approximation 
because it only effects the regions of poor sensitivity. 

The backgrounds for the point source analysis in Ref. [13] mostly come from atmospheric 
neutrinos for upgoing events, and from cosmic muons for downgoing events. For our anal- 
ysis, we assume that backgrounds can be suppressed to a sufficient level, which, for point 
sources, can be achieved by the angular resolution of about one degree for upgoing events. 
For example, if the photon counterpart is used to pre-select potential sources ("a priori 
source candidate search"), the number of background events is relatively low, see Table 3 in 



Ref. 13 . For downgoing events, the suppression of the cosmic-ray muon background is more 



difficult, which means that here the assumption of low backgrounds has to be interpreted 



"This is the factor weighting E dN/dE in the integrand of Eq. (12 1 if integration over Log(i?) is used. 
If the differential limit and flux are smooth enough on a logarithmic energy scale, the flux limit typically 
is below the differential limit since enough contribution to the integral is obtained. However, if spikes or 
sudden flux or differential limit changes are present, the flux limit may also exceed the differential limit 
locally. An example is the Glashow resonance process, which is sometimes taken out of the analysis for this 
reason. Note that in this study, Log(i?) — Log2o(i?) everywhere. 
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Figure 2: Neutrino effective area as a function of the neutrino energy for different declinations in 
IceCube-40. The thin sohd curves correspond to muon traclcs from j/^, the dashed curves to muon traclcs 
from Vr after r decay. The thick sohd curves show the neutrino effective area for either z/^ or if flavor 
equipartition between and v-y is assumed. Data taken from Fig. 8 in Ref. 13 and re-arranged. 
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with care and somewhat higher backgrounds are expected]^ In either case, we can reproduce 
the sensitivity for an E~'^ flux in Fig. 19 of Ref. [is] very well. We have also checked the 
flux of atmospheric neutrinos to be expected within one degree angular resolution, see, e.g., 



Ref. 41 . It turns out that especially in the energy range between 100 GeV and 10 TeV, 
where Earth attenuation effects are relatively small, the IceCube-40 point source sensitiv- 
ities for upgoing events are already touching the atmospheric neutrino background in the 
most conservative (systematics) case for the atmospheric neutrino flux. This means that in 
this energy range, we expect that the IceCube sensitivity will be background-limited very 
soon, whereas for the present analysis, neglecting the atmospheric neutrino background is 
a good approximation. 

We apply the above mechanism to obtain a limit for arbitrary fluxes to our model. We 
show in Fig. [s] the limits for selected u^^ spectra (including flavor mixing) for different 
declination bands in IceCube-40. We also show as thick lines the limits for an E~'^ flux 
(in the dominant energy rang^, and as thick curves the differential limits. The spectra 
correspond to the test points as marked in the figure, which correspond to the parameters 
R and B in Fig. [T] and a = 2. From the differential limits, we can read off that the optimum 
sensitivity for downgoing events is found at higher energies ~PeV, whereas upgoing events 



are best constrained at about 10 TeV 13 . The quasi-horizontal events exhibit a relatively 
broad differential limit. As far as the absolute sensitivity is concerned (see, e.g., thick solid 
lines for integrated limits), horizontal and upgoing events have similar sensitivities, but the 
downgoing events face lower statistics because of the necessary cuts to reduce the cosmic-ray 
muon background. As far as the limits for the individual spectra are concerned, we make a 
few observations. First of all, the best sensitivity, i.e., lowest normalization, for a particular 
spectrum is typically obtained if the spectral peak in E'^dN/ dE coincides with the minimum 
of the differential limit. The optimum sensitivity is then, to a first approximation, given 
by the maximum proton energy, which roughly determines the position of the peak and 
can be read off from Fig. [Tj For example, TP 11 is best constrained by upgoing events. 
However, especially if strong magnetic field effects are present and therefore the spectral 
shape becomes more complicated, such as for TP 2 for downgoing events, this rule does 
not necessarily apply anymore. In this case, there can be quite some impact of neutrinos 



from kaon decays, as we will discuss in Sec. 4.3, Note that, for a particular source, the 



declination is, of course, pre-determined. The way to interpret Fig. [3] therefore goes as 
follows: if a particular spectral shape for a source of declination 5 is described by a TP of 
this model, the sensitivity can be read off from the corresponding panel of Fig. |3} 

From Fig. [3| one can easily see that it is not trivial how to compare two different spectra 
with different spectral shapes. Consider, for example, the fluxes 2 and 11 in the upper left 
panel, both leading to the same event rate by definition. Which of the two neutrino sources 



^In fact, the cuts for the downgoing events seem to be chosen such that the backgrounds are roughly a 



factor of two higher for the downgoing events than the upgoing events, see examples in Table 3 in Ref. 13 



However, the absolute number of expected background events is nevertheless only a few, which means that 



the impact on the final sensitivity limit using the Feldman-Cousins approach 40 is relatively small on a 
logarithmic scale. 

^The choice of this energy range is somewhat arbitrary, as long as the main contribution to the integral, 
i.e., the minimum of the different limit is well contained in the integration range. Often the range where 
90% of the events are expected is shown. 
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Figure 3: Limits for selected t'^ spectra (including flavor mixing) for different declination bands in 
IceCube-40 (90% CL). The numbers of the individual curves correspond to the test points in Fig. [I] The 
thick lines show the limits for an flux (in the dominant energy range), and the thick curves the 

differential limits. Here the injection index a = 2 is chosen. 
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leading to these fluxes can be better constrained by the downgoing events? In order to 
quantify this aspect, it is useful to assign a single number to each spectrum which measures 
how much energy in neutrinos can be tested for a specific spectrum and event type. We 
choose the energy flux density 

as this quantity, which we show in units of erg cm ^ s ^ for point sources in order to distin- 
guish it from E'^dN/dE in units of GeV cm~^ s~^ (1 erg ~ 624 GeV). This quantity measures 
the total energy flux in neutrinos, and it is useful as a performance indicator measuring the 
efficiency of neutrino production in the source. 

In order to see that, consider the transformation of the injection spectrum of the neutrinos 
Q'^, (in units of GeV~^ cm~^ s^^) from a single source into a point source flux (in units of 
GeV~^ cm~^ s~^) at the detector dN/dE, which is (before flavor mixing and neglecting a 
possible Lorentz boost or beaming) given by (see, e.g., Ref. l33j) 



Here, V is the volume of the interaction region and dilz) is the luminosity distance. For 
the energy flux density, one has 



Audi 



where L^ = V j E'Q'JE' (15) 



is the "neutrino luminosity". Since the neutrinos originate mostly from pion decays and 
take a certain fraction of the pion energy (about 1/4 per produced neutrino for each charged 
pion), the neutrino luminosity is directly proportional to the (internal) luminosity of protons 
Lint (or the proton energy dissipated within a certain time frame AT) and the fraction of 
the proton energy going into pion production, commonly denoted by /^r- This quantity is 
a measure of the efficiency of pion production]^ Since a possibly emitted photon flux can 
be linked to Lint by energy equipartition arguments, one has oc x Lint oc f,, x L^, 
and is a measure for the pion production efficiency times luminosity of the source (if no 
photon counterpart is observed), or even the pion production efficiency itself (if a photon 
counterpart is observed). For example, for GRBs, this can be nicely seen in App. A of 



Ref. 42 , from which it is also clear that redshift and F cancel if the neutrino flux is 
related to a photon observation (the neutrinos and photons will, to a first approximation, 
experience the same Lorentz boosts, beamings, and redshifts). Therefore, the sensitivity to 
(j) for different spectra yielding the same number of events (as the fluxes do a the sensitivity 
limit) can be regarded as the prime performance indicator. Consider, for instance, two 
very similar sources producing neutrino spectra with slightly different shapes. Then the 
sensitivity to 0, and therefore to the pion production efficiency, may be very different, 
and it is fair to say that one source can be better constrained than the other. Consider, 
on the other hand, two sources producing very different neutrino spectra. In that case, 
the sensitivity to could be even similar, which means a similar f^^ x Lint is required for 

^It can be roughly approximated as /tt ~ 0.2R/Xpj for the A-resonance, where Ap^ = l/(n^cr) is the 
mean free path of the proton, n-, is the photon number density, and a is the interaction cross section. 
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detection. Our sensitivities to in the following sections have to be interpreted in that way. 
Finally, note again that our method does not predict the expected level of a neutrino flux, 
which means that in specific cases /tt x Lint might actually be much higher than in other 
ones for astrophysical reasons. Here, however, here we discuss only the interplay between 
spectral shape and detector response. 



4 Constraints from IceCube 



Here we present our main result, the sensitivity to the spectral shape of sources as a function 
of the Hillas parameters, and we discuss the dependence on the injection index a, the impact 
of kaon decays, and detection via muon tracks. 



4.1 Sensitivity on Hillas plane 



For a comprehensive parameter space scan, we compute for each set (i?, B) in Fig. [Tjthe z/^ 
spectrum including magnetic field and flavor effects. Then we normalize it with Eq. (12) 
for each event type, as it is illustrated in Fig. [3] for several examples, and we compute 
the energy flux density according to Eq. (13). The resulting sensitivity to (j) is shown in 
Fig. |4] as a function of R and B for a = 2 and the different declination bands in the 
different panels. The darkest regions mean highest sensitivities, as it is shown in the plot 
legends. For downgoing events, sensitivities as low as 10~^ ergcm"^ s~^ can be achieved, 
for upgoing events, sensitivities as low as 10~^''ergcm~^ s~^. The best parameter space 
coverage is actually obtained for near-horizontal upgoing events [5 = (0°, 30°)], which means 
that sources with relatively low x Lint can be detected there. 

As one result, the detector responds very well to the usual suspects, such as AGN cores and 
jets (TP 3 and 4), and to sources on galactic scales, such as TP 9 and 10. For instance, in 
the upper right panel (downgoing events), these are just in the optimal region. However, 
the best absolute sensitivities are obtained for TP 2 and 11 rather than TP 3 and 4, see 
upgoing events in lower row. The reason is the relatively low optimal energy for upgoing 
events, see Fig. |2} In addition, TP 5, 7, and 8 are not within the optimal sensitivity ranges, 
since these spectra peak at relatively high energies (see, e.g., Fig. [3] for TP 8). In fact, in 
our model, exactly these spectra stem from very high energy protons, i.e., these test points 
may be the best candidates to produce the highest-energetic cosmic ray protons (see regions 
marked "UHECR" in Fig. |4]). This just visualizes which might be an intrinsic feature of 
IceCube, which should be rather model-independent: the differential limits in Fig. [2] peak 
at relatively low energies, while the neutrino energies for interactions of 10^^ eV protons in 
the source may be expected at 10^^ to 10^° eV (10^° to 10^^ GeV). Therefore, IceCube may 
not be the best instrument to test the nature of the highest- energetic cosmic ray sources - 
although it may compensate for that by its size. Note that this argument does not depend 
on the composition of the UHECR, since it depends on the energy of the particles, not on 
the region of the Hillas parameter space. Moreover, there is a part of the parameter space, 
which is very difficult to test: especially sources with B > 10^ Gauss from downgoing events 
may be difficult to find in IceCube, which means that galactic sources with strong magnetic 
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Figure 4: Energy flux density [erg cm~^ s~^] sensitivity as a function of R and B for different declination 
bands in IceCube-40 (90% CL). The different contours (colors) correspond to the regions where a specific 
energy flux sensitivity can be exceeded. Here the injection index a = 2 is chosen. The dashed regions 
"UHECR" indicate where 10^" eV cosmic ray protons are expected to be produced in the model. 
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Figure 5: Flux limit for several different injection indices a = 1.5, 2, 2.5, and 3 for test point 4 (AGN 
jets) and two different declination bands (90% CL). The upper thick curves show the differential limits. 

fields are per se difficult to access, unless they are very luminous]^ Here DeepCore may 
improve the sensitivity substantially, since it has a lower threshold. 

4.2 Dependence on injection index 

So far, we have chosen a = 2 as the universal injection index for protons and electrons, 
as it may be roughly expected from Fermi shock acceleration. This, however, is not the 
spectral index of the neutrino spectrum. In our model, the target photons are produced 
by synchrotron radiation from the co-accelerated electrons, and the (synchtrotron) energy 
losses of the electrons are taken into account. As a result, the spectral index of the neutrino 
spectrum is approximately a/2 up to the maximal energy cutoff or the critical energy in 
Eq. ([5]), where the secondary spectra become loss-steepended by synchrotron losses. This 



spectral shape is similar, for instance, to neutrinos from blazar jets, see, e.g., Ref. 43 . We 
show a corresponding example for TP 4 in Fig. |5]for two different declination bands. In this 
figure, the sensitivities and differential limits are shown for several values of a. One can 
easily see that the spectral index of the neutrino spectrum between about 1 GeV and 1 PeV 
is roughly a/ 2. Since energy losses of the secondaries are small in this case, the spectral 
shape is relatively simple. For the case of the downgoing tracks (left panel), the peak of the 
spectrum coincides with the minimum of the differential limit. The energy flux sensitivity 
decreases as a increases, since more energy can be hidden in the neutrino spectrum at low 
energies. For the upgoing events (right panel), the peak of the spectrum is found at higher 
energies than the minimum of the differential limit. In this case, increasing a means that the 
neutrino spectrum can be better constrained and that the sensitivity, which is determined 
by the high energies, improves. 

It turns out that all test points in Fig. [l]can be separated in two categories: TP 1, 2, 11, 
and 12 always follow the trend in Fig. |5| left panel, since the neutrino energies are relatively 
low. The other test points exhibit a behavior similar to TP 2: the optimum sensitivity 

^Notc that the energy losses of the secondaries decrease the fraction of energy going into neutrinos even 
further. 
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Figure 6: Energy flux density sensitivity as a function of the injection index a for two selected test 
points (90% CL). The different curves represent the six different declination ranges used for IceCube-40, 
from downgoing [6 — (—90°,— 60°)] (dotted) to upgoing [6 = (60°, 90°)] (solid) with decreasing dash gaps. 
The gray-shaded regions mark the a-ranges which may be roughly plausible for Fermi shock acceleration. 

depends in tlie declination of the source. We show the energy flux sensitivity as a function 
of the injection index a for two examples (TP 2 and TP 4) in Fig. |6] The different curves 
represent the six different dechnation ranges used for IceCube-40, from downgoing [S = 
(—90°,— 60°)] (dotted) to upgoing [6 = (60°, 90°)] (sohd) with decreasing dash gaps. In 
the right panel, we flnd the functional dependence for TP 4, which we have qualitatively 
described above for upgoing and downgoing events. Note that the actual dependence on 
a is relatively flat for the downgoing and quasi-horizontal events, at least for reasonable 
injection indices (gray-shaded range), whereas for upgoing events the sensitivity signiflcantly 
improves with For TP 2 (left panel), low injection indices are always preferred, but the 
dependence on a is moderate. 

4.3 Impact of neutrinos from neutron and kaon decays 

Here we discuss two neutrino fluxes often not taken into account: For sources optically 
thin to neutrons, most of the neutrons which are produced in photohadronic interactions 
(see Eq. (|6])) will decay into electron antineutrinos either inside or outside the source. The 
protons produced in these decays may contribute to the cosmic ray flux. The neutrinos carry 
only a very small fraction of the neutron energy, which means that this flux is normally 
present at very low energies. However, for sources with very strong magnetic flelds, the 
neutrino flux from neutrons may actually be dominant, since the parents are electrically 
neutral. In addition, kaons can be produced by interactions such as Eq. ([7]). While the 
contribution of the neutrino flux from kaon decays is usually small, the higher value of the 

^At about a = 4, a minimum is found in these cases. At this value, the neutrino spectral index is 
roughly two, which means that the slope in the E'^dN/dE plot changes at around this value from increasing 
to decreasing. 
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Figure 7: Impact of neutrinos from kaon decays: Limits for several spectra (including flavor mixing) 
for two different selected declination bands in IceCube-40 (90% CL). The solid curves include neutrinos 
from kaon (and neutron) decays, the dashed curves represent neutrinos from the pion decay chain only. 
The numbers of the individual curves correspond to the test points in Fig. [l] The thick curves show the 
differential limits. Here the injection index a = 2 is chosen. 

critical energy Eq. ([s]) compared to muons and pions may lead to a dominant neutrino flux 
from kaon decays at high energies for strong enough magnetic fields. 

All these effects are fully taken into consideration in our calculations. Nevertheless, we 
show in Fig. [7] the impact of these two neutrino fluxes at two examples. In this figure, the 
limits for several spectra for two different selected declination bands are shown. The solid 
curves include neutrinos from kaon and neutron decays, the dashed curves show the limits 
if neutrinos from the pion decay chain only are considered. In the left panel, one can see a 
clear enhancement of the flux at high energies in all cases, which comes from the additional 
kaon decay component. For TP 2, the additional hump coincides with the differential limit 
minimum, which means that it contributes to the sensitivity (the solid curve is below the 
dashed curve for lower energies). The same applies to TP 1 and 12, although not so clearly 
visible. In the right panel, the contribution for TP 2 is negligible, since the solid and dashed 
curves coincide. There is, however, some contribution in the other two cases. For TP 1 one 
can also see the contribution from neutron decays dominating at low energies. However, 
while kaon decays help in parts of the parameter space, especially for downgoing events 
(where the differential minimum is at higher energies), the impact of neutron decays is, in 
general, small. 

4.4 Impact of muon tracks from 1/,-? 

In Fig. [2| we have compared the neutrino effective areas for muon tracks from and v^-. 
Assuming flavor equipartition between and z/,-, are there parts of the parameter space 
where the muon tracks from z/^ actually limits the sensitivity? In fact, we show in Fig. |8] 
several examples where we find some impact. In this figure, the thick solid curves show the 
differential limits including the Ur events (corresponding to the thick solid curves in Fig. [2]), 
the thick dashed curves the contributions from i/^ only. The thin solid curves show selected 
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Figure 8: Impact of v^. detection: Limits for several spectra (including flavor mixing) for two different 
selected declination bands in IceCube-40 (90% CL). The solid curves are based on the detection of muon 
tracks from both and Vr assuming flavor equipartition, the dashed curves are based on the detection of 
muon tracks from z/^ only. The numbers of the individual curves correspond to the test points in Fig. ^ 
The thick curves show the differential limits for both and Vr (solid) and i^^ only (dashed). Here the 
injection index a = 2 is chosen. 

spectra using tlie inclusive effective area, the thin dashed curves show the corresponding 
spectra for based events only. For the downgoing events (left panel), there is especially 
some impact at low energies, where the sensitivities can be significantly improved. For the 
upgoing events, the improved effective area for high energies has hardly any impact, at least 
for a = 2, for which the spectrum is parallel to the differential limit. Only for a < 2, some 
improvement may be expected. Re-drawing Fig. |4] including the z/^ events, the result mostly 
deviates a bit in the high energy region in the lower right. The effect at low energies is 
hardly visible, because the absolute sensitivities for TP 1 and 12 are worse than the shown 
contours. 



5 Complementarity to other experiments 

Here we qualitatively point out the complementarity among different data and different 
experiments, and we comment on the potential of future experiments. There are a number 



of high-energy neutrino data, such as from AMANDA 44 , 45 , IceCube-22 cascades [46] , 
Auger [20|, RICE |47], ANITA |48j, to name a few examples. Most of these data (except 
from AMANDA) have been applied to diffuse flux limits, which makes a comparison to 
the IceCube point source results, which are discussed in this study, per se difficult. While 
AMANDA has, in principle, a lower threshold than IceCube, it is not clear from the present 
literature if the neutrino effective areas are significantly better at low energies than for 
IceCube-40, since the larger volume may compensate for that. The RICE and ANITA 
experiments are based on radio detection initiated by cascades in the Antarctic ice. In 
these for cascades in IceCube, it is difficult to quantify the performance for a 

particular flavor. For example, in Ref. [46] (IceCube cascades), a flavor composition of 
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40% electron neutrinos, 45% tau neutrinos, and 15% muon neutrinos was given for an 
extragalactic test flux; see also discussion in the Appendix I in Ref. [47j (RICE). The only 
exception in this list is Auger: Earth-skimming v-r may produce tau leptons in the Earth, 
which may escape the Earth after energy losses, decay in the atmosphere, and produce 
an extensive quasi-horizontal, slightly upgoing air shower detectable by the Auger surface 
detector. Since this signal is practically flavor-clean, we illustrate our main points with 
2004-2008 data from the Auger experiment |20j, which appears to be simplest example. 
However, note that also the other experiments imply complementary information. 



5.1 Earth-skimming neutrinos in Auger 



The sensitivity for arbitrary fluxes for Auger can be very easily obtained from Eq. (12). We 



use the "conservative systematics" exposure from table III of Ref. 20 for our simulation. As 
the main differences to IceCube, diffuse flux limits are discussed, and the event sample is free 
of backgrounds from the atmospheric neutrino flux in that energy range. Since the viewing 
window constantly changes, there is a strong declination dependence of the sensitivity. 
However, for the sake of comparison, one may compute a "quasi-point source" flux by 
multiplying with the diffuse flux with the viewing window solid angle dVt ~ 0.6j^ Since this 
number is of order unity, we use the diffuse flux in the case of Auger directly. However, 
one should keep in mind this qualitative difference when one compares the numbers, which 
cannot be easily avoided]^ 

We show in Fig. [9] the limits for selected Vr spectra (including flavor mixing) for Auger 
2004-2008 data (90% CL). The numbers of the individual curves correspond to the test 
points in Fig. [l] The thick lines show the limit for an E'"^ flux (in the dominant energy 
range), and the thick curves the differential limit. Comparing to Fig. [sj one can easily see 
that especially the spectra peaking at very high energies can be very well constrained. Take, 
for instance, TP 8, and compare the result to Fig. |3j lower left panel. It is obvious that 
the normalization can in this case be better constrained by Auger, and, consequently, the 
energy flux density. We have also tested the impact of neutrinos from kaon decays here. 
Because the neutrinos from kaon decays show up at high energies and the Auger sensitivity 
(see differential limit) is dominant at about 10^ GeV, it turns out that the kaon component 
can be especially important in that case, see also Ref. |28|. One example, TP 3, is shown 
in Fig. [9} where the rightmost hump is the additional contribution from kaon decays which 
clearly limits the sensitivity. 
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Figure 9: Limits for selected spectra (including flavor mixing) for Auger 2004-2008 data (90% CL). 
The numbers of the individual curves correspond to the test points in Fig. [T] The thick line show the limit 
for an E^'^ flux (in the dominant energy range), and the thick curve the differential limit. Here the injection 
index a = 2 is chosen. 



5.2 Complementarity among different data and experiments 



In order to illustrate the complementarity of different data and experiments, we show in 
Fig. 10 the regions where the sensitivity exceeds 10~^ erg cm~^ s~^ [sr~^] for z/^ detection 
(left panel) and v.^ detection (right panel) for several selected data samples. In the case 
of IceCube, only muon tracks from or have been used, see neutrino effective areas in 
Fig. |2} In the case of Auger, the diffuse flux limit for Earth-skimming tau neutrinos has been 
used. Note that (at least if the r track cannot be separated) IceCube cannot distinguish the 
original flavor of the neutrino leading to a muon track, while Auger can. With this figure, 
we can qualitatively illustrate a number of points: 



Viewing window complementarity 

Different data samples correspond to different viewing angles. In the case of IceCube, 
the upgoing events test the northern sky, the downgoing events the southern sky. From 
the left panel, we clearly see the better sensitivity of, for instance, the quasi-horizontal 
upgoing events compared to the downgoing ones. For almost vertical downgoing events, the 
sensitivity is even poorer (c/.. Fig. |4]). New experiments built in the Northern hemisphere 
will complement this sensitivity, such as ANTARES or KM3NeT. 

^The solid angle d£t ^ 27r / dcos9 ^ 0.6, since 6 is integrated from 7r/2 to tt/2 + am with am — 0.1; see 
discussion after Eq. 4 in Ref. I20]. 

point source analysis at^Auger is presently being performed fio*, using 2004-2010 data. For a more 
detailed analysis than the one presented here, including the strong declination dependence, the exposure as 
a function of energy and declination would be needed, which is currently not publically available. However, 
from the comparison with Ref. [49] , our sensitivity with the 2004-2008 data roughly corresponds to a point 
source sensitivity at the declination 50° ^ l^| ^ 55° with the updated 2004-2010 data set (replace units by 
GeVcm^^s^^ on the vertical axis in Fig.^. 
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Figure 10: Regions where the sensitivity exceeds 10"^ erg cm^^ s^^ [si'^^] for i/^ detection (left panel) 
and v-r detection (right panel) for several selected data samples (90% CL, a = 2). The dashed regions 
"UHECR" indicate where 10^° eV cosmic ray protons are expected to be produced in the model. See main 
text for details. 



Energy range complementarity 

Different experiments or data will test different energy ranges. For example, from the right 
panel, IceCube can cover a wide region of the parameter space. However, for the high B 
(low energy) region, which is almost empty, experiments with a lower threshold will be 
needed. Especially the DeepCore part of IceCube is expected to extend into this region. 
In the lower right corner marked "UHECR", however, where the highest-energetic cosmic 
ray sources are suspected for protons, in fact Auger provides the best sensitivity, given 
our assumptions. Therefore, Auger and the radio-detection experiments may in fact be 



promising techniques to reveal the nature of these. Especially the Auger North project 50 



may achieve a substantially higher exposure in exactly that energy range 51 , as well as the 



JEM-EUSO project 52 . 



Flavor complementarity 



In Fig. 10, the left and right panel correspond to different neutrino flavors at the detector. 
If, for some reason, the equipartition between and is severely perturbed, it is first of all 
interesting that IceCube can, at least for upgoing or quasi-horizontal upgoing events, still 
test most of the parameter space relatively well (right panel). Of course, the Auger limit 
applies to only, but in the case of flavor equipartition it can be directly applied to z/^. 
It is now an interesting theoretical question when one may expect a strong deviation from 
the equipartition between z/^ and z/,-, which is a consequence of nearly maximal atmospheric 
mixing, at the detector. 
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In the Standard Model including massive neutrinos, flavor equipartition between these two 
fluxes relies on two necessary conditions: ^13 = and ^23 = 7'"/4, which guarantee flavor 
equipartition regardless of the initial flavor composition at the source. Note that exact tri- 
bimaximal mixings (sin^ ^^12 = 1/3) are not necessary. We have chosen these values in this 
study, which means that all of our results are exactly v^-Vt symmetric. However, there may 
be deviations from these values. Therefore, we have checked that for arbitrary initial flavor 
compositions (without Vr contamination at the source) the ratio between z/^ and Vr can, 
at most, vary between about 0.5 and 2 for the currently allowed 3o" ranges of the mixing 
parameters 



39 . This means that within the Standard Model, the and Vr fluxes at the 



detector have to be equal up to a factor of two, and any limit on (z/^) can be directly 
translated into a limit for v^. (z/^) within a factor of twoj^ On logarithmic scales, flavor 
equipartition is therefore guaranteed. 

This picture may change if physics beyond the Standard Model (BSM) is considered. First 
of all, strong perturbations of the flavor equipartition between and Vr at the detector 
may be unlikely if the physics BSM causes effects at the production or propagation of the 
neutrinos, since maximal atmospheric mixing will lead to equilibration of the two flavors 
again. Nevertheless, such exotic scenarios have been discussed in the literature, such as in 
the context CPT violation, see Ref. 53 for a discussion. Any new physics interaction at 
the detection (or, for upgoing events, in the Earth) may be most plausible reason to perturb 
this ratio. Examples are non-standard interactions at the detection, such as e^^, which are, 
however, limited to about 10% compared to the Standard Model [54] . Another possibility 
is a possible superluminal motion of muon neutrinos, OPERA has recently claimed at the 



6cr confidence level 55 , if a flavor dependent effect is at work which advances or delays one 
flavor compared to the others (see also Ref. [56]). In addition, it is not clear how the deep 
inelastic scattering cross sections evolve to high energies. Therefore, considering different 
flavor may still be considered complementary. 



6 Summary and conclusions 

The main motivation of this study has been the discussion of the interplay between spectral 
shape and detector response at neutrino telescopes beyond a simple assumption for 
the neutrino flux, from a particle physics perspective. Several effects often not taken into 
account have been included, such as magnetic field effects on the secondaries and flavor 
mixing. Especially, the impact of neutrinos from neutron and kaon decays has been dis- 
cussed, and the impact of r decays into muon tracks in the detector. As data samples, we 
have used 2008-2009 data from IceCube-40 for time-integrated point source searches, and 
2004-2008 data from Auger. As parameter space of interest, we have used the Hillas plot, 
described by the parameters R (size of the acceleration region) and B (magnetic field). For 
the description of the neutrino spectra, we have used a self-consistent model for the neu- 
trino production, in which protons interact with synchrotron photons from co-accelerated 
electrons (positrons) [21 1. While this model certainly does not apply to all spectra on the 

-•^^The deviations from flavor equipartition between and Vr can in fact be largest if the initial flavor 
composition is dominated by v^- 
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Hillas plot, it has been useful to illustrate some of the main qualitative points as a function 

of R and B. 

In order to compare the detector response to different spectral shapes, we have used the 
sensitivity to the energy flux density. The energy flux density has been introduced as 
a measure for the source luminosity L times pion production efficiency /tt (if no photon 
counterpart is observed), or even /^r directly (if the neutrino spectrum is compared to a 
photon observation). Here the pion production efficiency is the fraction of proton energy 
dumped into pion production. This means that the "sensitivity" in this work can be roughly 
interpreted as the sensitivity to f^^. Conversely, a source in a region of high sensitivity will 
be easier detectable, i.e., for lower luminosities and /^r, than a source in a region of low 
sensitivity. Note that we do not make any prediction for the expected level of the neutrino 
flux, as it can only come from astrophysical arguments. 

Since a particular source will be seen with a speciflc declination in IceCube, the flnal (shape) 
sensitivity will be declination-dependent. As one result, we have conflrmed that IceCube 
has excellent sensitivity to spectral shapes corresponding to AGN cores or jets for the model 
considered, which are among the usual suspects for signiflcant neutrino production, for all 
declinations. For example, they are found to be within the optimal region for downgo- 
ing tracks in IceCube. However, even better absolute sensitivities are obtained for quasi- 
horizonal and upgoing events for the spectral shapes from potential source classes such as 
white dwarfs, whereas the AGN blazar spectra may peak at a bit too high energies for these 
source declinations. The strong magnetic fleld region, for which the maximal proton energy 
is synchrotron hmited, cannot be easily tested by IceCube because of the relatively high 
threshold. Therefore, galactic sources with high magnetic flelds and signiflcant neutrino 
production are per se difficult to find. Data from DeepCore may signiflcantly improve this 
parameter region. Note that some optimization of the detector response may be performed 
for speciflc spectral shapes, which we have not taken into account. 

The parameter space region from which the highest-energetic cosmic ray protons may be 
expected, i.e., E 10^° eV, is in fact somewhat better covered by Auger, at least in princi- 
ple. This can be understood from generic arguments, recovered in the model used: 10^° eV 
protons colliding with much less energetic target photons will lead to neutrinos of about 
10^ GeV to 10^*^ GeV, which is just around the optimum of the differential limit of Auger. 
The conclusion holds for neutrino spectra signiflcantly harder than E^'^, which is a common 
expectation of the target photons come from synchrotron emission. Therefore, the common 
picture of an E~'^ neutrino flux can be clearly misleading, and the search for neutrinos from 
the sources of the highest-energetic cosmic rays may greatly beneflt from future experiments, 
such as Auger North or JEM-EUSO. This argument does not depend on the composition 
of the highest-energetic cosmic rays, since it is related to their energy, not their nature. 

We have also tested the dependence on other model parameters, such as the injection index 
a, which has some impact, especially for spectra extending to high energies for upgoing 
events in IceCube. Again, it is clear that general conclusions on astrophysical sources 
cannot be drawn independent of the neutrino spectral shape because of the interplay with 
the detector response. We have, for instance, shown several examples where the neutrino 
spectrum peaks off the energy range IceCube is most sensitive to, which means that rather 
large f^^ x L are required in order to be detectable. In turn, this means that a source 
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may not be found in spite of a correct astrophysical prediction of the flux, often based on 
energy equipartition arguments, simply because the flux shape does not match the detector 
response. 

As far as the particle physics is concerned, we have shown that neutrinos from kaon decays 
improve the sensitivities in cases where magnetic fields lead to an additional kaon decay 
hump in the spectrum and this hump coincides with the optimal sensitivity range. This 
effect can be especially prominent in Auger, since the kaon decay part shows up at the high 
energy end of the spectrum. In IceCube, the contribution from i/T--induced muon tracks 
(via the leptonic r decay channel) improves the sensitivities in parts of the parameter space 
if flavor equipartition between i/^ and at the detector is assumed, which means that it 
should be not a priori neglected. 

Finally, we have emphasized the complementarity among different event samples and ex- 
periments, such as with respect to viewing window, accessible energy ranges, and different 
measured flavors. For example, we have demonstrated that even if the equipartition between 
and Vr is strongly perturbed, IceCube can cover most of the discussed parameter space 
at least for upgoing events already from the i/^-induced contribution to the muon tracks. 
While such a perturbation can only be expected within the (neutrino) Standard Model up 
to a factor of two for the current mixing parameter uncertainties, new physics effects may 
be in charge of larger deviations. 

We conclude that the interplay between spectral shape and detector response is important 
for the detection of astrophysical neutrino sources. While the common assumption of an 
E^"^ neutrino flux is well motivated for GRBs, it is well known not to apply to different 
source classes, such as AGNs. Since the detector response strongly depends on the shape 
of the neutrino spectrum, the amount of energy which is needed to be dumped into pion 
production at the source to guarantee a neutrino detection may have to be actually higher 
than assumed in many astrophysical source models. 
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